Shot noise in chaotic systems: "classical" to quantum crossover 
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This paper is devoted to study of the classical-to-quantum crossover of the shot noise in chaotic 
systems. This crossover is determined by the ratio of the particle dwell time in the system, Td, to the 
characteristic time for diffraction tE — ^ ^ \ lnh\, where A is the Lyapunov exponent. The shot noise 
vanishes when tE 3> Td, while reaches a universal value in the opposite limit. Thus, the Lyapunov 
exponent of chaotic mesoscopic systems may be found by shot noise measurements. 
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The shot noise in the current flowing through a system 
is a random process associated with the nonequilibrium 
state into which the system is driven by the applied volt- 
age [|lj . The zero- frequency power spectrum of the noise 
is given by 



5': 



dt'[(I(t)I(t'))- (J) 2 ] =2Fe(I), 



(1) 



where — e is the charge of the electron, and (/) is the 
average current flowing through the system. In vacuum 
tubes, the shot noise emerges from the Poisson distribu- 
tion of the transmission of uncorrelated electrons. The 
result, in this case, is given by Schottky § formula (§), 
with F = 1. However, in more general cases, the charge 
carriers are correlated, for example due to Fermi statis- 
tics. These correlations suppress the noise, compared to 
Schottky result, by a factor known as the Fano factor, 
F. In diffusive wires with noninteracting electrons |}]]4|] 
F = 1/3, while in quantum dots F = 1/4. These 
suppression factors are universal in the sense that they 
are independent of the details of the systems. 

In a view of this universality, an intriguing feature of 
these results is the absence of an explicit fi dependence 
in F. This is despite the fact that the source of the shot 
noise is quantum mechanical uncertainty. Indeed, it has 
been shown by Beenakker and van Houten H that the 
shot noise vanishes in the classical limit. The purpose 
of this Letter is to characterize the quantum-classical 
crossover of the shot noise in noninteracting mesoscopic 
systems. It will be shown this crossover is governed by 
the ratio of two time scales: the average dwell time, Td, 
which is the typical time the electron stays in the system, 
and the Ehrenfest time, tE oc |log?i| which is the time 
scale where quantum effects set in. 

A qualitative understanding of the quantum-classical 
crossover of the shot noise for Fermions can be achieved 
by examining the general formula j^pl at zero frequency 




FIG. I. The classical-to-quantum crossover of the shot 
noise in chaotic dots. The solid line shows the Fano factor 
(^) as function of the ratio between the dwell time of the 
particle in the dot, Td, and the Ehrenfest time ([]). 

and temperature: S = Tr{T(l- T)} e 3 V/2Trh. Here V 
is the applied voltage, and T is the transmission matrix. 
For a single channel, the factor 7~(1 — T) represents the 
variance of the probability of a particle to pass through 
the channel. It vanishes both when this probability is 
unity (T = 1) or when it is zero (T = 0). Thus, the 
shot noise of open or closed channels is zero because the 
incident flux of Fermi particles, at zero temperature, is 
noiseless. The only possible source of the noise is the 
probabilistic nature of the transmission and reflection, 
which implies T / 1,0. 

Classical dynamics, however, is deterministic. A tra- 
jectory of an incident particle can either pass through 
the system (T = 1, if we loosely view the trajectory as 
a channel), or be reflected back (T = 0). Therefore, the 
shot noise vanishes in the classical limit. Yet, a quantum 
particle can switch between trajectories, an event which 
we term as diffraction ]8j . Consequently, the wave packet 
of the particle splits into two components: one passing 
through the system and another which is reflected back 
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(see inset of Fig. f ). The transmission probability is now 
different from one or zero, and the shot noise is finite. 
Thus the quantum-classical crossover of shot noise is de- 
termined by the probability of the particle to: (a) split 
its classical path by diffraction, and (b) that the result- 
ing two trajectories diverge during the dwell time of the 
particle in the system. 

Consider a chaotic system with smooth potential en- 
ergy characterized by a length scale a. If a is larger than 
the particle wavelength, \p, the typical time for the tra- 
jectories to diverge is the time which takes for a minimal 
wave packet to spread over a distance of order a. We 
picture the initial electron wave packet as a collection of 
trajectories separated by a distance Xf [see Ref. [|| for 
similar estimates for weak localization]. The separation 
among these trajectories after time t is dominated by the 
unstable nature of the classical dynamics. Namely, if p(0) 
is the initial distance between a pair of trajectories, then 
after time t the distance is of order p(t) ~ p(0)e At , where 
A is the Lyapunov exponent of the system. Thus, setting 
p{t E ) = a and p(0) = A^ one obtains: 



t E 
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Notice that the Ehrenfest time, t E oc | In 7a | , diverges 
logarithmically when h — > 0, while the average dwell time 
of the particle, Td, is essentially independent of Ti. Thus, 
in the classical limit, the Ehrenfest time is always larger 
than the average dwell time, t E 2> Td, and the shot noise 
is zero. The universal results for the shot noise men- 
tioned above, implicitly, assume the system to be in the 
opposite limit, 3> t E . In this regime the particle stays 
in the sample long enough to experience diffraction, and 
as a result the shot noise is finite, see inset of Fig. 1. 

So far, our discussion applies to any chaotic system. To 
make it more specific, consider the shot noise of quantum 
dots having a non-integrable shape. Let us assume the 
elastic relaxation time within the dot to be much shorter 
than the Ehrenfest and dwell times, and the contacts with 
the leads to have an equal large number of channels. Let 
the leads be in thermal equilibrium, and the bias between 
them be V . As we will show the zero-frequency correla- 
tor, at zero temperature, is given by Eq. (|l|) with the 
Fano factor: 
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where A2 is of order A. The behavior of this Fano factor 
is depicted in Fig. I. 

Before turning to the rigorous derivation of Eq. (^) , we 
discuss the physical meaning of this result. Factor F con- 
sists of the quantum value, 1/4, multiplied by the proba- 
bility of the particle to diffract during the dwell time, see 
also Ref. Indeed, for the chaotic dots the distribu- 
tion of the dwell time, t, is Poissonian, P{t) = e~ t ' Td It&, 



where Td is the average escape time. As we already ex- 
plained, only the trajectories with dwell time lareger than 
t E contribute to the noise. Thus, the relative number of 
such trajestories is e~ tE ^ Td which is the first term in the 
exponent of factor T in Eq. (Q). The second factor in the 
exponent take into account the deviations in the value 
of the Lyapunov exponent SX for different trajectories, 
(<5A 2 ) = A 2 /t. 

We turn, now, to the formal derivation of the above 
results. We choose to work in units where fi = k E = 1. 
The dependence on these constants will be restored in 
our final expressions. The derivation jl0| is based on the 
Keldysh formalism [fj"T| . 

The matrix Green function of our system 



<?(C;C) 



G R (C,C) 



<? A (C;C) 



C = (r,t) (4) 



(R, A, K denote retarded, advanced, and Keldysh com- 
ponents respectively) satisfies the Schrodinger equation: 



d_ 

dt 



H}1- a r A ■ J 



= i*(C-C). 



(5) 



where H is the Hamiltonian of the system, cr x A(() ■ J 
is a source term, a x is the first Pauli matrix in Keldysh 
space, and J is the current density operator: 



2m \dr 
The variational derivative 
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gives the current density correlation function. Here- 
inafter, £? j = AB + bA, for any operators A and B. 

Integrating Eq. (||) over the cross sections of the leads 
yields the total current correlation function (0). 

Our main purpose is to reduce the above equations 
to simpler ones which hold in the semiclassical limit, 
a 3> Xf- However, since the transition bewtween classi- 
cal trajectories plays a crucial role here, it is difficult to 
implement periodic orbit theory [l2| for this purpose, see 
e.g. O). The standard diagrammatic technique |l4| is 
suitable neither, because it is technically difficult to take 
the classical correlations into account. To circumvent 
this difficulty, we follow Ref. || and add a weak random 
potential, V(r), to the semicalassical potential U(r). The 
random potential, V(r), generates a small angle scatter- 
ing that models the diffraction; the total Hamiltonian of 
our system is H = H + V(r), where Ho = p 2 /2m + U (r) 
is the bare Hamiltonian of the system. To mimic the 
diffraction effects, the strength of V(r) is chosen such 
that the transport mean free time for the scattering on 
V(r) is given by 
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l/r tr = \X F /a. 



(7) 



The numerical coefficient in Eq. (Q) is not important for 
the calculation with logarthmical accuracy; value of r tr 
enters the final result only througth t F = (1/A) ln(Art r ), 
see Refs. |j 16 for further discussion. 

The construction of the semiclassical approximation of 
the Green function follows the usual steps |0,15[. First, 
the operators are Wigner transformed: 



g w (K,P;t,t') = J dre lPr g (r 
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Next, Q is averaged over the disorder in the self-consistent 
Born approximation, and the Green function is pro- 
jected on the energy shell. The projection is obtained 
by changing the coordinates from (R, P) to (£,x) where 
£ = 7Yo(R, P) — £f is the distance along the energy axis in 
phase space, and x = (R, n) are coordinate on the energy 
shell, n being a unit vector in the momentum direction. 
The outcome of the procedure is that the on-shell matrix 
Green function, 



ff(x;t,t') = — / "d£ £ w (£,x;i,t') 

■KV J 



(8) 



where v is the Weyl density of states, satisfies the 
Boltzmann-like equation: 



— — t 

dt + dt' + 
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g(x;t',t') = I[g}. 



(9) 



Here, the Liouville operator is given by 

t = dH 1 d__dH 1 d_ 

dV <9R <9R <9P' 1 ' 

while the source operator O is 

Og = iev F [A(x, t')ga x - A(x, t)a x g] , v F = v F n, 

and vf is the Fermi velocity of the electrons. Finally, the 
probabilistic effects in the problem are introduced by 
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J[ff] = -2^-[V»S*S-S*V£s], V n = nx — 

where r tr is given by (Q), and the convolution of two 
functions, say / and g, is defined as: 

(/*<?)(*,*') = ( dt f(t,t)g(t,t'). 



The homogeneous Eq. (^) is supplied with the constraint 

g*g = 5(t-H). (11) 

In order to find (|^), Eqs. ([)] - [ll]) should be solved in 
first order in A: g — go + g\ where g$, g\ are zeroth and 
first order in A respectively. The zeroth order term is 
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1 s(t-t>) r 



0: 



where the Keldysh component satisfies the equation 

ff( f(x;i,e) =0. (13a) 



Ot T tr 



The leads are assumed to be in thermal equilibrium, 
therefore, the boundary conditions are 



g (x — > ±oo; e) = 2tanh 



e ± eV/2 
2T 



(13b) 



where T is the temperature. 

Calculating the off-diagonal components of 



where g(e,u) = J dtdTg(t + ,t^)e teT+tuJ \ t± = t±r/2, we 
first solve the equation for the component gf ; 



-icj + C-\ V* 

Ttr 



gf = -2ieVFA u 



(14a) 



Then, we find the diagonal components from Eq. ( |11| ) as 
2g^ = —gff * gf , 2g^ = <j>f * gff , and substitute them 
into the equation for the component g^ . We thus obtain 



-iuj + C V* 

Ttr 



g^ = 2ievpA u 



(14b) 



-T^V n [ 5 ^(x;e+^). 9o K (x;e)V n5 f(x,^)] . 

Solving Eqs. ([lj), and substituting the results in Eq. (g), 
we find with the help of Eq. M) : 
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'-{{n a n' l3 (F 1 +F 2 ))), (15) 



where the LHS represents the uj component of the Fourier 
transform with respect to t — t', and ((• • •)) denotes an 
angular averaging with respect to n and n'. The entries 
in Eq. (p~5h are 



(16) 
Q(x ). 



Ft = X> w (x,x')Q(x')+X>-o ) (x',x)Q(x) 
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dxo 



O w (x,xo)D_ a ,(x' ) x ) 
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where Q.4 is the surface of unit sphere in d dimensions, 
2\>( x , x') is the classical propagator of the system, 
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P w (x,x') = O d( 5(x-x / ), (17) 



and Q characterizes the available energy space for the 
electron- hole pairs: 
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Q(w;x) 



de 



1 



g£(x;e+uj)g£ (x; e) 



K i 



(18) 



In equilibrium, g^ = 2tanh(e/2T), and consequently 



eq(<jj) = 2htU COth 



TlLU 

2T' 



(19) 



Thus, F2 vanishes, and Fi reproduces the fluctuation- 
dissipation theorem: i*i = 2Q eg ReP a ,(x, x'; w). 

Equations (jl^ - [l7]) describe the noise in any system 
with a semiclassical potential. To obtain more explicit 
results, we have to specify a model. We assume that the 
average dwell time, r^, is much larger than the classical 
ergodic time within the dot, and the contacts have the 
same number of channels. Moreover, we consider only 
noise at frequencies uTd <C 1 not to be concerned with 
effects of dynamical screening. 

The subsequent consideration follows the lines of the 
calculation of weak localization corrections ||. The re- 
sults for the factor ( ji~8| ) inside the dot (see Sec. IV of 
Ref. ||) take the form Q = Q eq + Q neq , where 

/r 

Q neq (tu, V) = —— ^2 lQe q (" ± eV) - Q eq (to)} , (20) 
± 

and r is given by (|). Within the same approximation, 
Q in the leads is given by the equilibrium value (19). 

The appearance of the factor W in © is not ac- 
cidental. In the absence of the small angle scattering, 
<7(^(x) can take only two possible values: either the equi- 
librium value of the left lead or that of the right lead. 
This is because a point x in phase space has a unique 
trajectory connecting it to incoming trajectories from the 
leads. Consequently, Eq. ([l|) implies that Q(x) = Q eq . 
However, diffraction, modeled by small angle scattering, 
allows for two distinct classical trajectories to reach the 
same point. At such a point, g$ (x) is a linear combina- 
tion of the equilibrium values of the leads, and therefore 
Q(x) assumes a nonequilibrium value. Thus, the factor 
yT 7 expresses the probability of two classical trajectories 
to become close to each other so that a transition from 
one to another due to diffraction is possible. 

Having the function Q neq in the form (^0|), the prob- 
lem of using Eqs. ( |l5|) and (|l6|) for the evaluation of 
the noise spectrum becomes equivalent to the evaluation 
of the weak localization correction (Sec. VI of Ref. Il) 
up to the replacement of the Cooperon with Eq. (|20|): 



C(l,l) 
trum: 



2~KvQ neq . We, thus, obtain for the noise spec- 



ked) = G 



Q eq {tu) + ~VF<3 



neq 



(21) 



where G is the conductance of the dot and Q eq is given by 
Eq. @. Substituting Eq. @ into Eq. @, we obtain 
the final result for the noise spectrum in the quantum 
dot with symmetric contacts: 



leg 



eV 

U) ± — ) - Q eq (w) 



where F is given by Eq. (g). In the limit T = 0, 
Q eq {u) = 2h\uj\ and S(ui) at ui = reduces to Eq. ([!]). 

To summarize, we constructed a theory for the 
quantum-to-classical crossover of the shot noise in meso- 
scopic system. A key ingredient of this crossover is the 
divergence of classical orbits in chaotic systems. This di- 
vergence is determined by the Lyapunov exponent. Thus, 
measurements of the of the shot noise in quantum dots as 
function of the dwell time can be used to determine the 
Lyapunov exponent of the underlying classical system. 
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